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1. Introduction 

It is of great challenge to understand a systematic pattern of heavy quarkonium mass spectra [|I|] 
from nonperturbative QCD. Lattice QCD simulations offer strong tools to gain an understanding of 
nonperturbative QCD, which enable us to compute the mass spectra of quarkonia directly by eval- 
uating appropriate correlation functions of operators on the lattice. However, due to a hierarchical 
mixture of ultraviolet and infrared physics in the quarkonium system, it is not straightforward to de- 
duce the underlying mechanism responsible for the mass spectra from the conventional simulation 
results. 

One way to solve the hierarchical problem in the quarkonium system is to employ effective 
field theories (EFTs). In this context, the use of nonrelativistic QCD (NRQCD) [Q, ||] and potential 
NRQCD (pNRQCD) ||, ^] have been proposed. The aim of EFTs is generally to have better 
control of a hierarchy of energy scales in an original theory for a high precision calculation. For 
instance, suppose there exist two typical energy scales which satisfy £"high 3> E\ ow in the original 
theory, then by integrating the scale above £"high> an EFT with the energy scale E\ ow can be derived. 
The high energy contribution is to be incorporated in the effective couplings of the interactions in 
the EFT, which are referred to as the matching coefficients, at any desired fixed order of perturba- 
tion theory in the small ratio -Eiow/^high "C 1. 

In the quarkonium system, it is considered that there exists a hierarchy of three types of energy 
scales, m S> mv ^> mv 2 , with a heavy quark mass m 3> Aqcd and a quark velocity v. NRQCD has 
been derived by integrating the energy scale above m in QCD, and pNRQCD by integrating fur- 
ther the energy scale above mv in NRQCD. Nonrelativistic nature of the quarkonium system then 
becomes manifest in these EFTs, and notably, a quantum mechanics-like hamiltonian emerges in 
pNRQCD. The matching coefficients in pNRQCD are dependent on the distance between a heavy 
quark and antiquark and eventually classified in powers of 1/m. Thus a set of these matching coef- 
ficients can be regarded as the heavy quark potential consisting of a static potential and relativistic 
corrections. Therefore, once these matching coefficients are determined, various properties of 
quarkonia, such as the mass spectra, can be investigated systematically by solving the Schrodinger 
equation. 

The nonperturbative matching is crucial for obtaining the heavy quark potential in pNRQCD 
as the energy scale of mv can be the same order of the magnitude as Aqcd- For this purpose, the 
present authors have been performing lattice QCD simulations and obtained accurate data so far up 
to the 0(1 /m 2 ) corrections [|7|, § g 0, [Hp. The aim of the present report is then to demonstrate 



a spectroscopy analysis of quakonia in pNRQCD with the lattice QCD results; the mass spectra of 
charmonium and bottomonium are computed by solving the Schrodingier equation, and the effect 
of relativistic corrections, in particular, the 0(1 /m) and the 0(\/m 2 ) spin-orbit corrections are 
examined systematically. The pattern of the mass spectra is found to be in fairly good agreement 
with experimental data, in which the 0(1 /m) correction gives an important contribution. 

2. pNRQCD and the heavy quark potential 

We begin by describing how the hierarchy of energy scales in the heavy quarkonium system, 
m S> mv S> mv , is controlled along the derivation of NRQCD and pNRQCD and how the heavy 



2 



Heavy quarkonium spectroscopy in pNRQCD with lattice QCD input 



Yoshiaki Koma 



quark potential is defined nonperturbatively. 

The first step is to integrate the energy scale above m in QCD, which leads to NRQCD. The 
meaning of the integration is that an interaction of heavy quarks via gluons whose momenta are 
larger than m is no more visible as a dynamical one and a creation or an annihilation of a heavy 
quark and antiquark pair cannot be seen within the resolution of the resulting obscured vacuum. 
Then, the heavy quark and antiquark are treated as independent external fields, respectively. 

Let \jf an d Xc be a quark and an antiquark annihilation operators, respectively, the NRQCD 
Lagrangian density with the leading interaction terms of the heavy quark and antiquark is written 
as §1 

D 2 D 4 1 

^nrqcd = Y\iD Q + ^ + + (V -> Xc) - 4 V M V 

+ / oB DE-ED a-(DxE-ExD), 

+Y r (c F g—+c D g — — — +ic s g — 2 )y + (Y^Xc) 

+^Vx^+^^f^+^ t r^7^+^ t r<T rZc t ra Zc , (2.1) 

m L m L m L m L 

where B and E represent the color-magnetic and color-electric fields, and Do and D are the time and 
spatial components of covariant derivative. The Pauli matrix is denoted as a, and T a (a = 1,...,8) 
are the SU(3) color generators. The factors cp, cp>, cs, d ss , d sv , d vs , d vv are the matching coefficients, 
which are to be determined so as to reproduce the same quantity, such as scattering amplitudes, 
both in QCD and NRQCD at any desired fixed order in perturbation theory of a s = g 2 / (4tt). This 
matching is always performed perturbatively as m S> Aqcd- In this way, the high energy momenta 
of gluons are integrated into the matching coefficients, while the remaining gluons in NRQCD 
carry nonrelativistic low energy momenta. For instance, if the matching scale is chosen just to 
be m, the matching coefficients of the bilinear term have the forms, cp = 1 + j= {Cp + Ca) + 0{a 2 ), 



cd = 1 + j^Ca + 0{a 2 ), cs = 2cp — 1 [ 12], where Cp = 4/3 and Ca = 3 are the eigenvalues of the 



quadratic Casimir operators of the fundamental and adjoint representation in SU(3), respectively, 
and those of the four-Fermi contact terms are d ss = ^na s + 0{a 2 ), d sv = —\n<x s + 0{a 2 ), d vs = 
— \na s + 0{a 2 ), d vv = \na s + 0{aj) Jl3|]. More general forms of the matching coefficients are 
provided in the original references [ |I2[ [Dp. 

The second step is to integrate the energy scale above mv in NRQCD, which leads to pNRQCD. 
A merit of further integration is to avoid a complication due to a mixture of the remaining two 
energy scales mv and mv 2 in NRQCD, which may affect a power counting of operators. The 
procedure is similar to the first step, i.e. one computes the same quantity in NRQCD and pNRQCD. 
However, in contrast to the first step, this matching must be performed nonperturbatively, since mv 
can be the same order of magnitude as Aqcd- The procedure proposed by Brambilla et al. [Q, [|] 
for the nonperturbative matching is called the quantum mechanical matching, which consists of the 
following steps. Firstly one writes the NRQCD hamiltonian with the \jm expansion, //nrqcd = 

+ 'h 1 ^^ H and then evaluates the expectation values of H^ 1 -^ with an eigenstate 

of H®' projected onto a color-singlet q-q state. Typically, the expectation values of are 
represented with color-field strength conelators on the q-q source. Then, the ground state energies 
of these expectation values are compared to the pNRQCD matching coefficients, which are also 
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classified in powers of 1 /m, 

Vrqcd = V (0) + ~V {1) + + ■■■ , (2.2) 

where V^ 1 - ' are local in time but nonlocal in space (dependent on the q-q distance r), hence they 
can be regarded as the components of the heavy quark potential. \A°) is expressed with the Wilson 
loop, and thus, related to the static potential. W'- 1 ) are clearly due to finiteness of the quark 
mass, which are expressed with the color-field strength correlators and identified as the relativistic 
corrections to the static potential. So far these expressions are provided up to of 0(l/m 2 ) H|. 



The second term in Eq. ( |2.2[ ), corresponding to the 0(1 /m) correction, is expressed as 



V^(r) = - dtt((E(0)-E(0))). (2.3) 



This term is typical in QCD due to the self interaction of gluons (three-gluon vertex). The double 
bracket {(&)) indicates that an expectation value of an operator G on a q-q source normalized 
by the expectation value of the q-q source itself. The third term contains the spin-dependent and 
spin-independent corrections. The spin-orbit correction in V^ 2 \ which is relevant in the following 
analysis, is expressed as [|]] 

y%\r) = ( g^T + y M Jdtt((Bj(0)E k (0))) +e ijk r i f dt t ((B\0)E k (r))))\l ■ s . (2.4) 

= V[(r) = V&t) 



The arguments of the color-field strength tensors in Eqs. ( |2.3| ) and ( |2.4[ ), such as or r, mean that 
a color-field strength tensor is attached to either quark or antiquark. The variable t is a relative 
temporal distance between the two color-field strength tensors, which is to be integrated. The 
matching coefficients V^-°'(r) are then computed by lattice QCD simulations. 

3. Brief summary of the lattice QCD results 

The present authors have been studying the heavy quark potential in pNRQCD by using lattice 
QCD simulations within the quenched approximation [Q, ||, ||, [l(], The numerical procedure 
has been to use the Polyakov loop correlation function for the q-q source, and compute the field 
strength correlators by employing the multilevel algorithm, and evaluate them with transfer matrix 
theory. So far the potential with the 0(1/ m ) relativistic corrections have been computed accurately 
up to the distance r ~ 1 fm. 

Fig. [| summarizes a part of the lattice QCD results which are relevant in the following spec- 
trum analysis, the static potential, the 0(1/ m) correction, and the 0(1/ m 2 ) spin-orbit correction. 
The Sommer scale ro = 0.50 fm has been used to fix the lattice spacing for each j3 value. The func- 
tional forms of the potentials from short to long distances have been determined by the following 
fitting functions, 

V(0)( r) = _^ + ar + c (0) ) ( 3. 1} 
r 

V (i) (r ) = _^ + a (i)l nr + c (D ) (3.2) 
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Figure 1: Summary of the lattice QCD results used in the present analysis, the static potential (upper 
left), the 0(1 jm) correction (upper right), and the 0(1 /m 2 ) spin-orbit corrections V[ (lower left) and 
V£ (lower right). Solid lines are the fitting curves. 
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The fitting parameters have been found to be 



a = 0.297 , a = 1 .06 GeV/fm , a (1) = 0. 142 GeV z , e = 0.2 . 



(3.3) 



(3.4) 



As can be seen in Fig. [j], the fitting curves nicely describe the behavior of the lattice data. Note 
that the leading order perturbation theory results in V^°\r) = — " with a = CfOC s , V^'(r) = ~^r, 
and V[(r) = 0, V 2 '(r) = a/r 2 . The short distance behavior of v(°), and V 2 ' may follow these 
functions, although the coupling a is different from the bare one. It can be observed that both 
V[ and V 2 ' contain long-ranged nonperturbative tails, which are related to the string tension in 
through the Gromes relation dV^ jdr = V{ - V{ @. 



4. Heavy quarkonium mass spectra in pNRQCD with the lattice QCD input 



The Schrodinger equation 

,2 



'2m + ^+V( >V {0) =£ ( V 0) 
, m J 



(4.1) 
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Figure 2: The mass spectra of charmonium for various quantum numbers, 2s+l Lj; a set of three red lines 
in each state from left to right are the spectra of E^°\ + AZi' 1 ', and E^ +AE^ + AE^\ respectively, 
while the black lines are the experimental data. The energy levels are normalized at the spin-averaged IS 
state. 




Figure 3: The same figure as Fig. 2, but for bottomonium. 



is then solved and the energy E^ and the wave function l// ' are computed. Using the wave function 
the corrections to Z?w are evaluated in the first order perturbation theory, AZJW = (y/"(°)|y 
for i > 1. In this way, the heavy quarkonia of various quantum numbers are investigated systemat- 
ically. The matching coefficients in NRQCD are dependent on the matching scale between QCD 
and NRQCD, and thus different for the charm and bottom quark sector, but they are simply as- 
sumed here to be the leading ones, cf = cs = 1, which correspond to taking the matching scale 
to be infinity. This issue will be discussed in the forthcoming paper. It should be emphasized 
that there is no free parameter except for the quark masses and a constant shift of the potential. 
The quark masses are chosen to be m c = 1.27 GeV (charm) and nib = 4.20 GeV (bottom) in this 
analysis. 

Figs. |2] and || show the mass spectra of charmonium and bottomonium, respectively, for various 
quantum numbers classified by 2S+x Lj. The energy levels are normalized at the spin-averaged IS 
state. A fairly good agreement can be observed between the computed pattern of the mass spectra 
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and the experimental data, especially for the states below the DD threshold for charmonium and the 
BB threshold for bottomonium, in which the effect of the 0(1 jm) correction is remarkable. As is 
clear from the I s operator, the 0(1 /m 2 ) spin-orbit correction affects the splitting among the levels 
of 3 P/=o,i,2 states, the fine splitting, where ( 3 Pj\l-s\ 3 Pj) = —2,-1,1 for J = 0,1,2, respectively. 
The splitting is less noticeable than that caused by the 0(1/ m) correction as the matching scale 
effect of QCD and NRQCD is not yet properly taken into account. The experimental level order of 
3 Pj states exhibits AE( 3 Pi — 3 Pq) > AE( 3 P2 — 3 Pi), but if only the spin-orbit interaction contributes, 
this will be opposite regardless the functional form of the spin-orbit correction due to the factor of 
( 3 Pj\l -s\ 3 Pj) . It would be crucial to take into account the spin-tensor contribution. 

5. Summary 

We have studied the charmonium and bottomonium mass spectra in potential nonrelativistic 
QCD (pNRQCD) with the lattice QCD results of the heavy quark potential. In pNRQCD, the 
heavy quark potential consists of a static potential and relativistic corrections classified in powers 
of l/m, and in the present investigation, we have examined the effect of the 0(1/ m) correction 
and the 0(\/m 2 ) spin-orbit correction systematically. We have found that the pattern of the mass 
spectra is in fairly good agreement with the experimental data, in which the 0(1/ m) correction 
gives an important contribution. This effect has not been taken into account in the conventional 
phenomenological models. It is certainly interesting to see all the effect up to the 0(1/ m 2 ) cor- 
rection in pNRQCD on the spectra, which contain the velocity-dependent potentials as well as the 
other types of the spin-dependent corrections. 
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